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Abstract 

We introduce a microscopic model for the dynamics of the order book 
to study how the lack of liquidity influences price fluctuations. We use the 
average density of the stored orders (granularity g) as a proxy for liquidity. 
This leads to a Price Impact Surface which depends on both volume lj 
and g. The dependence on the volume (averaged over the granularity) 
of the Price Impact Surface is found to be a concave power law function 
{(l>{u),g)}g ~ UJ^ with S « 0.59. Instead the dependence on the granularity 
is (j){u},g\uj) ~ g" with a —1, showing a divergence of price fluctuations 
in the limit g — + 0. Moreover, even in intermediate situations of flnite 
liquidity, this effect can be very large and it is a natural candidate for 
understanding the origin of large price fluctuations. 



1 Introduction 

One of the basic assumptions of the Arbitrage Pricing Theory, which is a stan- 
dard tool of Financial Engineering, is the perfect liquidity of a market [U E] ■ 
This assumption implies that a market is always able to absorb incoming or- 
ders without producing a significant movement in the price. However in the 
last ten years a great amount of studies demonstrates that finite liquidity ef- 
fects play a crucial role in order to understand the features of a modern market 
[H El [71 [HI E [ini ■ From these studies liquidity appears more important in fixing 
the price response of a market than the order volume or the flow of information 
which are canonically assumed in Economics as the main ingredients to explain 
price movements [TI] [T^ [T3] . 

The problem of finite liquidity have been highlighted by detailed studies of 
real order books [HI [TSl [HI . In order to explain the empirical features of or- 
der books, several analytical and numerical models have been introduced [17j . 
A very important class of models for the order book are the so called zero- 
intelligence models which are characterized by random order deposition mech- 
anisms [TFIITO]. 

In this paper we introduce and discuss a zero-intelligence model of order book 
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to investigate the relation between finite liquidity and price fluctuations. Zero- 
intelligence model are in general simple and workable but still able to reproduce 
the main properties of real order books. 

In order to quantify the liquidity of a given configuration of the order book we 
measure it via the average density of the stored orders in the book which wc call 
granularity. The analysis of the system response to a granularity fluctuation can 
be properly addressed in the framework of a really elementary dynamics for order 
placement which differs from the coarse-grained approach used in [THl [501 HI] ■ 
We find that the impact of an order, which depends both on its volume and on 
the granularity of the configuration of the book, is a concave function of the 
volume. This reproduces the empirical observations and represents an impor- 
tant test for the model itself. Our main result is that the price response to an 
incoming order is inversely proportional to the granularity. This means that 
granularity and the associated finite liquidity represent fundamental elements 
in the price response, which becomes divergent in the limit of vanishing gran- 
ularity. From these results it is clear that a realistic modelling of a financial 
market requires the inclusion of finite liquidity. 
The paper is organized in the following way: 

in section 2 we describe how an order book works and we introduce the problem 
of finite liquidity; 

in section 3 we define our model and point out similarities and differences with 
respect to other order book models; 

in section 4 we give the criteria to set the model parameters according to em- 
pirical data; 

in section 5 we introduce granularity as an operational measure of liquidity; 
in section 6 we introduce the Market Impact Surface and we study its depen- 
dence on both the volume and the granularity; 

in section 7 we draw the conclusions and outline the possible extensions of our 
model. 

2 Order book and the role of liquidity 

The elementary mechanism of price formation is a double auction where traders 
submit orders to buy or sell. The complete list of all orders is called the 'order 
book'. In Fig. [T]we give a schematic representation of an order book. There 
are two classes of orders: market orders and limit orders. Market orders corre- 
spond to the intention of immediately purchase or sale at the best price (quote) 
available at that time. The limit ones instead are not immediately executed 
since they are offers to buy or sell at a certain quote which is not necessary the 
best one. If we consider a sell limit order this means that its quote is higher 
than (or equal to) the best bid b{t) which is the order of buying with the high- 
est price. On the other hand a buy limit order implies that its price is lower 
than (or equal to) the best ask a{t) which is the order of selling with the lowest 
price. The non-zero difference between a(t) and b{t) is defined as the spread 
s{t) = a{t) — hit). The prices of placement of orders (called 'quotes') are not 
continuous but quantized in unit of ticks whose size is an important parameter 
of an order book. 

The price of a stock can be conventionally defined as the mid-price p{t) — 
\a{t) + b{t)]/2 and it can change only if a limit order falls inside the spread or if 
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Figure 1: Schematic representation of the order book dynamics. There are two 
classes of orders: the market orders and the hmit orders. The market ones are 
orders of purchase or sale at the best available quote. On the other hand the 
limit orders are not immediately executed since they are placed at a quote which 
is less favorable than the best quote. The volume lu of the orders is an integer 
multiple of a certain amount of shares. The price of a stock can be defined as 
the mid-price p{t) = [a{t) + b{t)]/2. 



a market order matches all the orders placed at the best quote. It is clear that 
the specific configuration of an order book it is a very important aspect for price 
movements. A thick book full of orders can absorb the arrival of new orders 
without giving rise to large jumps of the price. On the contrary if the book 
is sparse, even a small incoming order could trigger a large price variation. In 
order to clarify this problem we can consider an order book in the configuration 
of panel a in Fig. [2l This situation corresponds to a very liquid market in 
which the order book can absorb several orders without large price variations 
even if they are relatively large. This regime corresponds to the assumptions of 
the standard financial picture where an order can be immediately executed, its 
impact is marginal and the market is nearly efficient. 

Instead if a liquidity crisis occurs, the configuration of the order book changes 
dramatically and the situation is like the one represented in panel b in Fig. (5] 
The orders stored are few, the average distance between them is large and the 
fiow of market orders cannot be easily and immediately absorbed by the sys- 
tem and by consequence even a small order can produce a large price variation. 
Hence a lack of liquidity can dramatically amplify the response of the system. 
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Figure 2: A very liquid order book vs a liquidity crisis. This figure illustrates 
how the degree of liquidity of a market plays a key role in determining the 
system response to the volume of an incoming order. The order book is very 
liquid (panel a) when a great amount of orders is stored in each side of the 
order book and almost all quotes behind the best one are occupied. In such 
a situation a market order produces a small perturbation of the system and 
then a small price adjustment. On the contrary when a liquidity crisis occurs 
(panel b), the order book is characterized by few orders stored and by a large 
average gap between adjacent orders. In this case even a market order with 
a small volume can produce a dramatic price fluctuation of several ticks. In 
our discussion high and low liquidity situations are symmetric and refer to the 
structure of the order book. 



3 Definition of the model 

The purpose of this paper is to study the price response function in presence of 
liquidity crisis and to appropriately address this problem we are going to intro- 
duce a model with a suitable microscopic dynamics. Our model can be defined 
as a 'zero-intelligence' model as the ones introduced in (THl ^U\- It should be 
noted that the specific questions that we consider cannot be addressed by the 
dynamics of Refs. [HIED]. 

In fact in [121 US the deposition is modelled in term of a flow of orders driven 
by a Poisson process. The orders arrive at each quote at a certain fixed rate 
so that the deposition interval and the number of orders per unit of time are 
infinite but the orders per unit price interval are finite. The motivation of the 
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authors of [THl HO] is that in such a way the order book can be never depleted 
and some finite size effects, due to a finite deposition interval, are avoided. 
On the other hand in our model it is preferred an elementary microscopic mech- 
anism for order deposition. 

The difference between these two mechanisms for order deposition is the defi- 
nition of the time unit. In the present paper the time unit corresponds to the 
time (that is not fixed) between two following orders. Instead in Refs. [THl HH] 
the time unit corresponds to a physical and fixed amount of time (for example 
At equal to some minutes) during which several actions are made by investors. 
This permits a coarse-grained description in terms of average quantities, such 
as order flow. 

For the specific implementation of our model we start by considering three mech- 
anisms of deposition which will be compared with the empirical data. 
We want to point out that only one of the three mechanisms gives rise to a real- 
istic dynamics. Certain elements are however common to the three mechanisms 
and we start by discussing them. 

In our set-up, as in a real order book, an order can be placed at a quote that 
is an integer multiple (positive or negative) of the tick size Aq. In this paper 
we do not investigate the effects of this parameter; therefore we set it equal to 
I. For the sake of simplicity we assume that all orders (limit and market) have 
the same volume lu = 1. The fact that a market order has the same volume of a 
limit order is coherent with the observation that in more than 99% of the cases 
a market order does not exceed the best limit order available . 
In this model there is one mechanism for order creation, the limit orders de- 
position, and two for order annihilation, the arrival of market orders and the 
cancellation process of limit orders. The cancellation of a limit order occurs 
when the order has been stored in the book for r time steps without being 
executed, where r is a parameter to be fixecfl. The balancing of the three mech- 
anisms of creation/annihilation fixes the properties of the steady state that is 
reached very quickly (about 5 x 10'^ steps if t ranges between 100 and 1000). 
Two limit orders can be also removed from the order book when the spread 
s{t) is zero, that is when an incoming limit order is placed at the opposite best 
quote. However, these events are very unlikely (< 1% of the number of the 
market orders) and their effect is negligible. 

We define the price as the mid-price p{t) — [a{t) + b{t)]/2 where a{t) and b{t) 
are the quotes of the best ask and of the best bid. 

In our model we completely neglect the daily closures of the market and there- 
fore the strong price variations during the night. The typical length of a run 
is 10^ time steps (10^ if a larger sample is needed) which corresponds to a real 
sample of about 5 years, estimating about 10'^ operations per day. 

3.1 Order deposition 

Now we arc going to discuss three possible cases for order deposition which we 
will test with respect to empirical data in order to select the most appropriate 
mechanism. Due to the symmetry of the order book, the probability that an 
order is a sell or a buy one is always the same. 

^We know that this assumption, reasonable for the present study, is far to be reahstic from 
the moment that analysis of real order books have shown that the lifetime of a limit order 
increases monotonically with its distance from the mid-price |15l I21| 
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3.1.1 Case 1 



Once the nature (buy or sell) is determined, the order can be a market one 
with probability tt and a limit one with probability 1 — tt, with n < 0.5. Limit 
orders will be placed in the interval [b{t) , b{t) + L] or [a{t) — L, a{t)] depending on 
their nature (sell or buy respectively). The probability distribution within these 
intervals is considered uniform and L is a free parameter to be fixed. Hence the 
deposition interval is independent on the spread size. 

3.1.2 Case 2 

The order can be a market one with probability tt and a limit one with proba- 
bility 1 — TT with TT < 0.5 as in case 1. A limit order will be placed in the interval 
[b{t) + 1, b{t) + k s(t)] or [a{t) — k s{t), a{t) — 1] with k > 1 with a uniform distri- 
bution, being s{t) the spread. In this case the length of the deposition interval 
depends on the spread size through a sort of auto-regressive process. This im- 
plies that the probability that a limit order produces a price variation (i.e. an 
order is placed inside the spread) is time-independent and approximately equal 
to k~^. We anticipate that this second case will turn out to be the more realistic 
one with respect to the empirical data. 

3.1.3 Case 3 

The third case is inspired by the mechanism proposed in |21j . First a number 
^ in the interval [—L, L] is extracted. If ^ < s(t) the order is a limit order 
otherwise it is a market one. If the order is a sell limit order its quote is a(i) — ^. 
On the contrary if the order is a buy one its quote is b{t) + ^. In such a way the 
probability of being a market or a limit order is not set a priori. 

We have previously mentioned that our choice for order deposition produces 
a finite number of orders stored in the order book and a finite depth. This 
causes a finite size effect: the order book is subject to complete depletion with a 
cut off in those quantities linked with the book depth. However the probability 
of a complete depletion goes very quickly to when the cancellation parameter 
r is increased (this probability is less than 10^^ for r > 200). 
Moreover we note that the uniform deposition mechanisms used in this paper 
cannot reproduce in a realistic way all the empirical features of the order book; 
nevertheless this assumption has been made with the aim to define a minimal 
model still able to reproduce the main properties of the order book. Besides, 
the main goal of this work is the evaluation the effects of the granularity in price 
fluctuations, which depends on the static realization of the profile. 

4 Calibration of parameters and preliminary re- 
sults 

In this section we check the properties of the models we have introduced with 
respect to the empirical data. This analysis will also allow us to decide how to 
calibrate the parameters. 

The empirical results of Ref. [7] indicate that a realistic value for tt is approxi- 
mately 1/3. In case 3 the probability of being a market or a limit order is not 
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fixed a priori. Nevertheless in this case we find that the fraction of market orders 
tends to a reahstic value (about 0.32) when r grows, and becomes substantially 
independent on r when the probability of emptying the order book turns to be 
neghgible (i.e. t > 200). This feature does not seem to depend crucially on the 
parameter L. 

We now have to fix the remaining parameters. Two quantities that can be cal- 
ibrated on empirical data are the average spread and the average number of 
stored orders n. The parameters of the model which rule these quantities are 
r, k (case 2) and L (case 1 and 3). Since it is not always possible to a find a 
configuration of all parameters able to reproduces these two quantities in a real- 
istic way, we give priority to those sets of parameters that produce a reasonable 
average number of orders. 

All three cases exhibit a realistic accumulation of orders (i.e. 50 — 100 orders 
per side) for r ranging from 200 to 500. This clarifies the role of t in driving 
the average number of orders in the steady state. In case 1 we choose L = 200 
and in case Z L — 100, in such a way the order book has approximately the 
same depth in both cases. For these two cases we applied the previous criterion 
of priority in order to have a realistic average number of orders. Case 2 appears 
to be very interesting because it exists a range of parameters for which we can 
reproduce both a realistic average spread and average number of orders. For 
instance for r = 400, fc = 4 we obtain an average number of orders per side of 
about 60 and (s) « 2.8 ticks. In Table[l]we report a summary of the results for 
the tested parameters. We have only reported some selected results and a more 
detailed and systematic analysis will be given in a future survey |27j . 
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300 < r < 750 


200 < T < 500 
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/ 


3 or 4 
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General 
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remarks 


(s) > 20 


(s) = 2.8 


(s) > 20 






for A: = 4, T = 400 





Table 1: Summary of the setting of the parameters. Only case 2 is able to 
reproduce reasonable values for the quantities of interest. 



It is worth noticing that the auto-regressive case 2 is also interesting with respect 
to the Stylized Facts. In fact only in this case we can observe some volatility 
clustering in the price increments (returns). All cases exhibit a return probabil- 
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ity density function which decades as a power law (see Fig. [3]) but only in the 
second case the exponent 7 of the tail is compatible with the empirical observa- 
tions (—5 < 7 < —2). For instance in case 2 we find 7 = —4.1 for t = 200 and 
7 = -5.4 for T = 400. 



5 Operational estimator of the granularity 

Given a certain configuration of the order book we want to define a way to 
measure the liquidity of it. In fact our objective is a detailed study of the 
relation between finite liquidity and price fluctuations. 

We would like to describe the effect that the order book may not be compact 
but characterized by orders separated by voids. In order to perform this analysis 
we are going to introduce a function g which we call granularity. We propose a 
definition for granularity which is related to the inverse of the size of the average 
void between two adjacent quotes (gaps). At each time t the spread s(t) sets a 
characteristic length which can be used as unit of measurement for these gaps. 
However we are also interested in measuring the granularity of the order book 
in a region far from the best quotes since the two sides of the order book usually 
have a depth much larger than s{t). Therefore we define a partitioning of size 
s{t) of one side of the order book and we perform the following average 



, Nit) 



where N(t) is the number of intervals of length s{t) defined by the partition 
(we require N{t) > 2 in order to calculate g{t)) and ns(t),j is the volume of 
orders in an interval of length s{t). From the definition of N{t) we observe that 
N{t)s{t) = L is approximately the range of the order book and consequently 
the granularity g results equal to the average linear density of orders 

{git)) = J (2) 

where = X)^^*'' ''T-s{t).j- For the sake of simplicity we now drop the temporal 
dependence of g and the brackets of the average. We note that the definition of 
granularity given in Eq. Q does not measure the real average size of the voids 
because more than one order can be stored at the same quote. Therefore g~^ can 
be seen as the equivalent gap between two adjacent quotes in an hypothetical 
system where each quote can store only one order. We also observe that in the 
limit of a continuous order book, the price variation Ap produced by an order 
of volume lo would be 

f Ap 

p{p)dp = UJ (3) 



where p{p) is the density of the stored orders. If we approximate p{p) with its 
average value g we find the following scaling relation for g and lu 

Ap=-. (4) 
9 

In this framework a liquidity crisis occurs when 5 — > 0, on the contrary the 
market is very liquid if g ^ 1. In fact when g is very large a great amount of 
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Figure 3: Fat Tails. We show the probabihty density function of returns for the 
three mechanisms for order deposition and for different values of t. Because of 
the symmetry of the model the left tail is identical to the right one. We see 
that all three mechanisms give rise to a power law behavior but only the case 
2 produces realistic values for the tail exponent 7: 7 = —3.1 for t — 200 and 
7 = —5.4 for T = 400. In fact for cases 1 and 3 the exponent 7 is larger than 
—2 which is not realistic. Moreover we note that the case 1 and 3 are very 
similar with respect to the tail exponent and to the cut off around 100 due to 
the deposition mechanism. We conclude that these two mechanisms produce 
the same results except for slight differences. Conversely the auto-regressive 
deposition mechanism of the case 2 appears to belong to a completely different 
class. This fact is also supported by the fact that the case 2 reproduces a certain 
volatility clustering differently from case 1 and 3. 
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orders can be executed without producing a significant price variation. 
The granularity defined by Eq. ^ has the dimension of a price~^. In order to 
obtain an absolute parameter for liquidity, we define the dimensionless liquidity 
a.s g — gAq. As we have seen in section [3] at this stage we are not interested in 
the eiTect of the tick size Aq on the model so we have set it equal to 1. Con- 
sequently the dimensional and the dimensionless liquidity numerically coincide. 
Furthermore we want to point out that the definition of g given by Eq. is 
not directly related to a specific disposition in the order book but it gives only 
an average information about the void size. 

In Fig. m we show the histogram of g for the model labelled as 'case 2'. The 
probability to observe a liquidity crisis, as expected, vanishes when r — > oo 
because of the large amount of orders stored in the order book. However, this 
probability cannot be neglected for values of r ranging from 50 to 400 that we 
have previously recognized as the realistic ones. 

It is important to note that in this paper we give a different definition of gran- 
ularity with respect to the definition of Ref. [20]. In fact in that paper the 
authors define two dimensionless parameters to describe the discreteness of the 
system: the first, named granularity, is the dimensionless order size, the second 
one is the dimensionless tick size. In a certain sense both aspects are included in 
our definition of granularity and combined in a single parameter. These different 
choices are motivated by different purposes, while the authors of pO] gives the 
priority to the parameters that describe the order flows, differently we are going 
to focus our attention on how the discreteness influences the response function 
of the system. 



6 Price Impact Surface 

In this section we restrict our analysis to the 'case 2' mechanism for order 
deposition and all the results reported here correspond to r = 400 and k — 4. In 
order to analyze the average response of the system to an external perturbation 
(i.e. the volume of an incoming market order) after a certain time lag At, we 
introduce the Price Impact Function 77 deflned as 

Tjiuj, At) = 2E[Ap{At)\uj] . (5) 

One can consider two interesting limits of Eq. ([5]) , the asymptotic limit {At —^ 
00) and the instantaneous limit {At — > 0). In this paper we restrict our in- 
vestigation to the second case because the persistent effects on the market 
cannot be fully explained by a zero-intelligence model which neglects all the 
time-correlated structures of the order deposition (see[Hl [HI EHl [Ml HH])- 
Now we are going to consider also the dependence of </> on the granularity g. In 
this way we have to deal with a Price Impact Surface rather than a Price Im- 
pact Function and we can analyze how the price response of the system depends 
both on the order volume and on the granularity g. Hence we define the Price 
Impact Surface as 

(j){u,g)^2E[Ap\iu,g]. (6) 

Before investigating the dependence of the price response surface on g we should 
verify the well-known empirical property that {4>{uJ-, g))g (i.e. the standard Price 
Impact Function defined as in Eq. ([5])) is a concave function with respect to the 
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Figure 4: Probability density function of the granularity g for the case 2 for 
different values of r. Obviously large values of r allow for a significant accu- 
mulation of orders at every quote (i.e. g > 1), despite this fact the probability 
of g < 1 and consequently of large price variation is not negligible for r « 400 
which has been previously recognized as a suitable choice of this parameter. 



order volume ld. This is a manifestation of the fact that the order book is far 
from being in a linear regime. While the agreement on the concavity is almost 
universal, different functional form of the Impact Function have been given: the 
authors of [TS] propose ~ Inoj, instead other authors propose (/>(w) ~ lo^ 
with exponent 5 ^ 0.5 (^) and 5 « 0.3 ([3 EI])- In Fig. [Hwe find that the 
surface averaged on g (black crosses) is concave and follows a power law 

(0(tj,5))g - oj^ with 5 « 0.59. (7) 

6.1 Price Impact Surface in the direction of uj 

Now we consider the Price Impact Surface as a function of lo for fixed values of 
g. We see from Fig. [5] that its behavior is a power law where the exponent (3 
and the amplitude depend on g 

<P{u;,g\g)^Lo^. (8) 

This function is highly concave for small value of g and [3 grows for increasing 
values of g. However we notice that the dependence of the amplitude on g is far 
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Figure 5: Market Impact Surface in the direction of volumes w for case 2 and 
r = 400. We can see that the function is concave as the one observed for 
real data. Its behavior is approximately a power law uj^ . The dashed lines 
correspond to the results of the fitting procedure and they are shifted for clarity 
purposes. The price response for large values of g cannot be larger than the one 
for small values of g. The black crosses represent the impact function averaged 
on g and correspond approximately to an exponent of 0.59. 



stronger than the one of (3. This implies, as expected, that the Impact Surface 
produces smaller variation of price in the limit of large even if the exponent 
f3 is larger. 

It is worth noticing that the observed scaling behavior for lo is different from 
the one predicted by the zero order approximation made in Eq. |4l This is due 
to the peculiar shape of the order book profile whose maximum is peaked far 
from the best quote. 

The results of Fig. [5] are qualitatively similar to those of Fig. 6 of Ref. [5D] even 
if the different definitions of granularity do not permit a detailed comparison of 
the two models. 

6.2 Price Impact Surface in the direction of g 

Now we analyze the price variation as a function of the granularity g for a fixed 
volume oj of the incoming market order. We find that the g-dependence of the 



12 



surface is well represented by a power law 

(j){Lj,g\uj) ^ g" with a « —1. (9) 

In this case the exponent is negative and it appears to be nearly independent 
on oj (see Fig. [5]). In fact we find that its value is —1 (within the error bars) 
for all values of to. Consequently when the granularity vanishes the response 
of the system diverges. We can now quantitatively analyze how this limit is 
approached for finite (but large) values of g. For example if we consider an 
order of size w = 10, the average price variation induced when g w 0.5 is about 
ten times greater than the one observed when g 1. 

Therefore we have found that the amplification induced by the discreteness of 
the system is approximately proportional to the equivalent gap measured by g 
as predicted by the zero order approximation of section [5l However we want to 
stress that for small values of g some deviations from the g^^ scaling behavior 
seem to be observed. This fact suggests that, when g vanishes, the deviation of 
the density p(j)) from the average density measured by 5, is significant. Further 
investigation of this aspect will be considered in future works. 
The fact that the exponent of the power law in the direction of g does not seem 
to exhibit a significant dependence on uj may suggest a possible factorization of 
the price impact surface, as we are going to verify in the next section. 

6.3 Factorization of the Price Surface Impact 

In Ref. 15J the authors observe that the empirical impact function admits the 
following factorization 

7jiLo,At)^ fiioMAt). (10) 

Similarly we verify whether the Impact Surface Function could be expressed in 
the following way 

^u;,g) = ^{u;)ij{g) (11) 

To check if the factorization of Eq. (fTT|) is correct we have to verify if it exists a 
function of uj such that the impact surface divided by this function turns out to 
be independent on uj. This is equivalent to say that we would like to observe a 
complete collapse into a unique curve by rescaling with this function the surface 
(j). The simplest choice for this scaling function is the average impact function 
{(j){u),g))g ^ cj°-^^. We report the rescaled functions for different values of uj in 
Fig. [7] and we observe a nearly perfect collapse for w < 30 within error bars. 
In principle we cannot obtain a perfect factorization because of the dependence 
of the exponent (3 on g. Nevertheless we observe that this dependence is very 
weak, especially for g > 1 when (3 ~ 0.60. In fact in Fig. [7] one can appreciate 
some discrepancies only for g < \. We can conclude that the factorization 
proposed in Eq. (fTT|) approximately holds at least for small values of lo. We 
also expect that tf{Lo, At 0) can be factorized with respect to its variables, 
instead the situation may be very complicated for the function V'(ffj ~* 0) 
since g depend on time. We will consider these points in future investigations. 

7 Conclusions 

In this paper we have introduced a model of order book to study the generalized 
Price Impact Function (Price Impact Surface) and its dependence not only on 
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Figure 6: Market Impact Surface in the direction of the hquidity g for case 2 
and T = 400. In this direction the surface can be fitted with a power law t;" 
with a negative exponent. This implies that when g tends to (i.e. when a 
liquidity crisis occurs) the price response can be very large even if the external 
perturbation (i.e. the volume uj of an incoming order) is very small. The 
exponent a is approximately independent on oj and its value is compatible with 
— 1 within error bars. 



order volume but also on granularity. We find that the granularity operates 
as a strong amplifier of price variations when a liquidity crisis takes place. In 
particular the dependence of the Price Impact Surface on the granularity is a 
power law with exponent nearly equal to —1. This result implies that the system 
response to an incoming order diverges in the limit of vanishing granularity. It 
would be interesting to compare the Price Impact Surface found in the model 
with the empirical one. 

Furthermore agent-based models for financial markets usually do not take into 
account the problem of finite liquidity. In fact a common way to model the 
price movements in an agent based model is through the Walras' mechanism in 
which the price adjustments are proportional to the excess demand i.e. Ap = 
XED. The coefficient x ^-nd the excess demand ED are generally assumed to 
be independent on granularity. If, in first approximation, one interprets the 
Market Impact Function as the quantity the price movements are proportional 
to in the Walras' mechanism, then the previous coefficient x must depends on 
granularity. Moreover we have observed that the Price Impact Surface can be 
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Figure 7: Factorization of the Price Impact Surface. In order to verify if the 
Price Impact Surface admits the factorization 4>{uj, g) — ip{LL!)ip{g), we rescale the 
surface in the direction of g (i.e. 4>{uj^ g\uj) ) with a function of uj. For the sake of 
simphcity we choose the Impact Function averaged on g that we have previously 
found to be proportional to w'^-^^. We obtain a very good collapse of the rescaled 
functions within error bars. This confirms that the Price Impact Surface can be 
factorized at least for small values of In the insert we show a magnification 
of the rescaled functions to highlight their statistical compatibility. 
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factorized and we can try to identify the excess demand with tp{w) and x with 
ipig)- Consequently we aim at introducing this dependence on g in the workable 
agent based model we have introduced in previous papers [531 [Ml [13 [M] ■ In 
this framework we expect that even small unbalances of the market can produce 
large price adjustments if a liquidity crisis is present. Therefore the role of the 
amplification introduced by ip{g) could be one of the explanation of the breaking 
of the cause-eff'ect relation and so that small perturbations can produce large 
fluctuations. 
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